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ABSTRACT 

We consider dust drift under the influence of stellar radiation pressure during the 
pressure-driven expansion of an H ii region using the chemo-dynamical model MAR¬ 
ION. Dust size distribution is represented by four dust types: conventional polycyclic 
aromatic hydrocarbons (PAHs), very small grains (VSGs), big grains (BGs) and also 
intermediate-sized grains (ISGs), which are larger than VSGs and smaller than BGs. 
The dust is assumed to move at terminal velocity determined locally from the balance 
between the radiation pressure and gas drag. As Coulomb drag is an important con¬ 
tribution to the overall gas drag, we evaluate a grain charge evolution within the H ii 
region for each dust type. BGs are effectively swept out of the Hii region. The spatial 
distribution of ISGs within the Hii region has a double peak structure, with a smaller 
inner peak and a higher outer peak. PAHs and VSGs are mostly coupled to the gas. 
The mean charge of PAHs is close to zero, so they can become neutral from time to 
time because of charge fluctuations. These periods of neutrality occur often enough 
to cause the removal of PAHs from the very Interior of the Hii region. For VSGs, the 
effect of charge fluctuations is less pronounced but still significant. We conclude that 
accounting for charge dispersion is necessary to describe the dynamics of small grains. 

Key -words: dust, extinction - ISM: bubbles - Infrared: ISM 


1 INTRODUCTION 


Theoretical and observational studies of dust evolution 
and survival in Hll regions have a long history. The 
first evidence for dust within an Hll region was found 
from observations of scattere d light in the Orion Neb¬ 
ula do ’Dell & C. Hubba^ 1 19651 ) and several other H ll re¬ 
gions I O’Dell. Hubbard fc Peimbertlll96^ . These first ob¬ 
servations hinted that dust content in Hll regions de¬ 
creases towards ionizing stars. Development of infrared 
(IR) detectors has allowed the observation of not only 
light scattering b y dust in Hll reg i ons but also proper 
dust emission (e.g..lNev fc AllerJll96d : IStein fc Gillettiri969l : 

' 31 


iHaroer fc Low! 1 197^ . These observations have generally 

confirmed the presence of inner cavities with l ow dust den¬ 
sity i n compact and extended H 1 1 regions fe.g..lGillett et ^ 


19751: A itken, Griffiths & Jonei 
Chini. Krueael fc Wargauil98 
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Ii 


1977: Nakano et al. 


Aannestad fc Emery 


■ 19831: 

200 ill . 


On the other hand, evidence has also been pres ented in 
favour of dust overabund ance inside H ll regions dPanaeial 
1 1974 iTenorio-Tagld fl^l . 

The interest in the dust in H ll regions was renewed by 
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observations of Spitzer Space Telescop^ Herschel Space Ob- 
servator'^ and WISE. One of the great Spitzer results was 
a discovery of numerous IR emission bubbles in our Galaxy. 
IChurchwell et aP (l2006l l argued that these bubbles are pri¬ 
marily formed by hot young stars in massive star forming 
regions. Detailed study of the bubbles observed by Spitzer 
showed that 86 per cent of thes e objects indeed enclose H ii 
regions (iDeharveng et al1l20ldl . 

T he IR bubbles look like partial or closed rings on 8 pm 
maps llChurchwell et al]l2006l : ISimpson et al.llioid l. A sim¬ 
ilar ring-like morphology is s een in longer wavelen gth im¬ 
ages obtained with Herschel jAnderson et aLll2012ll . in the 
sense that in the far-IR the bubbles look like rings coinciding 
with 8 pm emission ring s or residing somewh at further away 
from the ionizing source JPaladini et al.l2012l l . An IR bubble 
morphology at 24 pm is significantly different. Emission at 
this wavelength is observed not only towards the outer ring 


^ Spitzer Space Telescope is operated by the Jet Propulsion Lab¬ 
oratory, California Institute of Technology, for the National Aero¬ 
nautics and Space Administration. 

^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 


© 2015 RAS 













































2 V.V. Akimkin, M.S. Kirsanova, Ya.N. Pavlyuchenkov, D.S. Wiehe 


but also towards the interior of an H ll region dWatson et alJ 
I 2 OO 8 I : [Peharveng et al.ll201Clf ). It must be noted that such a 
morphology was actually predicted bv I WrightI l| 1971^ 1 . 

The different appearance of H ll regions at different IR 
wavelengths must be somehow related to various factors 
(stellar wind, magnetic field, ultraviolet emission), which 
leave their imprints on dust properties in the vicinity of hot 
massive stars. To decipher this information, a detailed model 
is needed that takes into account not only relevant physical 
processes, but also a multitude of dust properties. Emission 
at 8 /rm is generally assumed to belong to PAH molecules ex¬ 
cited by ultraviolet (UV) photons, while emission at 24 pm 
is presumably generated either by stochastically heated very 
small grains (VSGs, a « 20 — 50 A) or by hot sub-micron 
grains (big grains, BGs). It is natural to attempt to relate 
the different spatial distributions of emission at near-, mid-, 
and far-IR wavelengths with spatial distributions of grains 
having various sizes and/or chemical compositions. 


While numerous models of dusty H ll regions, both 
static and dynamic, are presented in the literature, they of¬ 
ten either assume that dust is dynamically coupled to the 
gas or consider dust motion relative to the gas for a single 
dust component, being primarily directed towards estimat¬ 
ing t he ability of the dust to pa ss momentum to the gas 
(e.g., iKrumholz fc Matznerj l2009fl or to absorb ionizing ra¬ 


diation within the region (e.g., Petrosian. Silk fc Fieldl[l972l : 
[Arthur et al. 1 12003) . The first hydrod ynamic model of this 
kind was presented bv iMathewsI (Il967l l. who showed numer¬ 
ically that a central ‘hole’ in the Rosette Nebula can be 
explained by the radiation pressure on the dust, provided 
the dust is able to survive in the interior of an Hll region. 
He noted the importance of taking the grain charge into 
consideration, as a charged grain is more tightly bound to 
the gas due to Coulomb drag and is also protected against 
sputtering. Dust was assumed to be dynamically coupled to 
the gas. 


The relative motion of gas and dust u nder the influence 
of rad iation pressure was considered by iGail fc Sedlmavil 
lll979ll . These authors studied the formation of inner cavi¬ 
ties via the expansion of dusty cocoons surrounding embed¬ 
ded young massive stars, with the radiation pressure being a 
main driver for this process. Their model al lowed for the gas- 
dust d rift, treated as a diffusion process. iGail fc Sedlmavn 
1 I 1979 I) also stressed the importance of grain charging, point¬ 
ing out that the relative gas-dust motion is fastest in regions 
where the absolute dust charge is small, so that Coulomb 
drag is less effective. Variations of dust charge along the H ll 
region radius may lead to dust piling up at a certain range of 
distances from an ionizing star (or stars). As charging prop¬ 
erties depend on the grain size, density enhancements for 
big grains and small grains were shown to occur in different 
regions. 

A detailed analysis of radia tion pressure in a static 
dusty H II region was presented bv lDrairi^ (l2011ah . who con¬ 
structed a family of solutions for various H ll region proper¬ 
ties and estimated parameters of an inner cavity produced 
by radiation pressure when the dust is perfectly coupled to 
the gas. He also considered dust drift and concluded that 
dust grains can attain quite significant drift velocities (of 
the order of 100 km s”*^ and greater) in the inner part of 
the region, but this drift only affects a minor fraction of the 


entire dust content of the region. No distinction between 
various dust species was made in that study. 

Spatially resolved IR observations of Hll regions draw 
a complicated picture that apparently cannot be explained 
by the action of a single dominant mechanism on some 
generalized dust population. An ideal model would in¬ 
clude gas and dust motion, radiation pressure, stellar 
winds, grain charging, multiple dust populations, etc. We 
initiated our study of dust evolut i on in Hll regions in 
IPavlvuchenkov. Kirsanova fc Wiebd ll2013ll . It was shown 
that images of the RGW120 bubble at 8, 24 and 100 fim 
can be reproduced simultaneously if PAHs are partially de¬ 
stroyed inside the Hii region, but the abundance of VSGs 
and BGs stays close to its mean interstellar value. Dust emis¬ 
sion at 24 /rm in that model is mainly produced by stochas- 
tically heated VSGs as in the general interstellar medium 
llGompiegne et al.l l201ll ~l. while 100 fim emission is gener¬ 
ated by cold BGs in the dense shell surrounding the Hii 

region. _ 

_ The model from [ Pavlyuchenkov. Kirsanova fc Wiebd 

(l2013ll is based on the simplifying assumption that the dust 
is dynamically coupled to the gas and the only process that 
may alter the dust properties is the destruction of PAHs. 
In the present work, we study another, less explored aspect 
of dust evolution in an Hll region. Namely, we neglect dust 
destruction for the time being, but relax an assumption of 
tight coupling between the dust and gas and study how the 
differential dust drift under the action of radiation pressure 
and grain charging influence the spatial distribution of dust 
grains of different sizes and types within an H ll region. 


2 MODELLING 


The presented model is based on the chemo-dynamical code 
MARION, wh ich is being developed a s an extension of the 
Zeus-2D code (IStone fc Normanlll992fl . The MARION code 
is designed to simulate a pressure-driven expansion of an 
H II region around a young massive star. It is assumed that 
a star is embedded in a molecular cloud with a uniform gas 
distribution, characterized by number density n. At time t = 
0, the star starts ionizing and heating the surrounding gas. 
The size of the H ll region grows because of the gas pressure 
difference between hot ioniz ed and cold molecular gas (e.g. 
iMathew'illQGSl : [Spitzeilll978l l. A shock front develops ahead 
of the ionization front and propagates into the molecular 
cloud, b eing preceded by fronts of H 2 and CO dissociation 
(see e .g. iTielens fc Hollenbachlll 985l: iHosokawa fc Inutsukal 
I 2 OO 6 I : [Kirsanova. Wiebe fc Sobolev I^od l . This shock front 
shovels molecular gas into a dense envelope around the H ll 
region. 

The MARION code is desc rib ed in de- 
tail in Kirsanova. Wiebe fc Sobole v ([ 2 OO 9 II and 
[Pavlyuchenkov. Kirsanova fc Wiebd HioiSi l. Here we 
recapitulate its main features. The expansion of an Hll 
region is modelled in a one dimensional approximation, 
taking into account a detailed thermal balance of the gas 
and dust along with the related chemistry. In this work 
we do not consider the chemical evolution of the molecular 
gas, but we still have to account for the simple chemistry 
of H 2 , CO, H 2 O and OH molecules, which is needed to 
model heating and cooling processes in a photon dominated 
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region. Cooling in lines of neutral and ionized O, S and Si, 
which is important for the thermal balance of an H ll region, 
is also included. A summary of the thermal processes is 
given in Table [T] 

Photoreaction rates are calculated by in¬ 
tegration of corresponding c ross-sections, taken 
from the Leiden database dvan Dishoed3 1 19881 : 
[van Dishoeck, Jonkheid fc van HemertI 2006h . with the 


local emission spectrum used at each point of a computa¬ 
tion grid. Parameters o f other reactions are taken from th e 
UMIST 95 database llMillar. Farauhar fc Willacvl Il997ll . 
Initial physical parameters and non-zero abundances of 
chemical species are listed in Table [2] A full list of chemical 
species includes H, H+, Ha, CO, C, C+, O, 0+, 0++, S, 
S^, Si, Si”*", OH, H 2 O and He. Surface species are not 
taken into account to save computational time. Initial abun¬ 
dances in Table are based on the interstella r medium-like 
set o f abundances from the CLOUDY model llFerland et al.l 
l2013h. The spectru m of ionizing stellar emission is taken 
from iKurucd lll979h . 

For this work, we further upgraded the MARION code 
and added a treatment for charged grain motion relative to 
the gas. Two primary forces acting on a dust grain are radia¬ 
tion pressure and aerodynamic drag. The radiation pressure 
force is given by 


/- 




( 1 ) 


where c is the speed of light, F(v) is the radiation flux com¬ 
puted in the model of the H ll region, C^^iy) is the radiation 
pressure cross-section for the grain type g: 


Crpil') = + (1 - (cOs6l®))C'®ca(i^)- (2) 

Absorption and scattering cross-sections for carbonaceous 
and silicate dust grains, and Cl^a, are calculated us¬ 
ing Mie theory. PAH ab sorption cross-sections are adopted 
from iDraine fc Lil (l2007l L 

The grain charge is determined by the balance between 
electron a nd ion accreti on on to the grain and photoelectric 
emission (iTielend [200il (see Appendix El for details). To 
save computational time, the grain charge is pre-calculated 
on a grid of electron number densities, gas temperatures and 
distances from the star (or radiation field strength). 

A grain of radius Ogr and charge Zg^, moving with veloc¬ 
ity V relative to the gas, undergoes braking due to collisions 
with gas particles. Neglec ting surface processes, t he result¬ 
ing aerodynamic force is toraine fc Salpet^l 19791 ): 

Fdrag — 27mg,.A:T’gas ^ ^ [Go(^i)~b 

-I- Zi)G2{si)\ I, (3) 


where 

Si = 

G'o(si) « 8si/(3-y7r)y^l -|- 97rs?/64, 

G2(si) ~ Si/(3v^/4-l-sf), (4) 

/(ttgr X /bT’gas)^ 

A = 3/{2agre\(j)\)^/kT^Z/(^DD). 


Table 1. Heating and cooling processes. 


Process 


Reference 


Photoionization of H 
Photoionization of C 
Photoelectric heating 
H 2 formation 
H 2 dissociation 
H 2 pumping 


Heating 

Tielens (2005) 
Jonkhei^e^ah 

jJonkhei^^^aL 

Jonkhei^e^aL 

jJonkhei^^^aL 

Jonkheid et ah 


f2004^ 

f2004'l 

(2004) 

f2004'l 

[2004) 


Recombination of H"*" 

Lyman a emission 

Hydrogen free-free transitions 

[OI]A = 63.2^m 

[OI]A = 6300. SA 

[OII]A = 3728.8A 

[CII]A = 157.7Mm 

[SII]A = 6730.8A 

[OIII]A = 5 OO 6 . 9 A 

[OHI]A = 51.8//m 

[SIII]A = 953 O. 9 A 

[SIIIJA = 187.0/rm 

[Sill] A = 35.0/rm 

CO rot&vib transitions 

H 2 rot&vib transitions 

H 2 O rot&vib transitions 

OH rot&vib transitions 


Cooling 


Tielens 

^ielens 

Tielens 

^raine 

Tielens 


120051 
120051 
120051 
2011bl 
120051 
Hollenbacl^_&_^cKee 

Tielens 120051 
120051 
120051 
120051 
120051 
120051 

_ 2011bl 

Hollenbacl^_^^cKee 

Neufeld&Kauftn^ 

Hollenbacl^_&_^cKee 

Hollenbach fc McKee 


Tielens 

Tielens 

Tielens 

Tielens 

Tielens 

Draine 



Dust-gas interaction 


Tielens 


Here Tgas is the gas temperature, e is an elementary charge, 
k is the Boltzmann constant, rie is the electron number den¬ 
sity. Summation over i involves number densities m and 
masses rrii of all neutral and charged gas components. We 
consider H, H 2 , H'*’, He and electron number densities. Their 
radial profiles are taken from the results of chemo-dynamical 
simulation at a considered time moment. 

The gas velocity Ugas is given by the hydrodynamic part 
of the MARION code. We assume that a grain always moves 
with the terminal velocity Vgr, which is given by the grain 
equation of motion, provided that the grain acceleration is 
zero. In other words, we assume that as the grain moves 
through the H ll region, its velocity quickly adapts to a local 
environment. In this approximation 

Fdrag(Ugr, Zg-^^V = Vgas ^gr) “t“ Frp(Ugr) = 0. (5) 

Solving this non-linear equation relative to Vgr for grains of 
different types, one can simulate the motion of various dust 
grains in the Hll region. We have tested this assumption 
for integrating dust equations of motion in a static Hll re¬ 
gion pre-computed with the CLOUDY code llFerland et al.l 
l2013h and found that the time needed to reach the terminal 
velocity almost never exceeds 200 yr. 

In the numerical model we define Ugas and Vgr at cell 
borders and use them to calculate the corresponding fluxes 
across borders. Equation is solved with the bisection 
method for every cell. As soon as all Vgr are found, we use 
them to calculate dust fluxes fgr in the s ame way as in the 
trans port step of the hydrodynamic solver (IStone fc Normanl 
I 1992 II . The calculated fluxes are used to update the dust 
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Table 2. Initial conditions for modelling of an expanding Hii 
region. 


Quantity 

Value 

Tgas 

10 K 

Tdust 

From balance equation 

n(H + H2) 

3 X 10® cm“® 

a:(S) 

3.24 X 10“® 

a;(Si) 

3.16 X 10"® 

:r(H2) 

0.5 

a;(CO) 

2.51 X lO”'* 

a;(H20) 

6.8 X 10“® 

:r(He) 

9.8 X 10-2 


density according to the equation of mass conservation: 

1 a , 2 , S r, ICi\ 

^ + = ( 6 ) 

In this formulation, dust is dynymically coupled to the gas 

when Ugr — Ugas- 

Dust drift is considered on top of the hydrodynamic 
model, in the sense that the motion of grains is modelled 
simultaneously with the Hll region expansion, but without 
feedback. Specifically, to model dust drift and to calculate 
dust densities we use the following approximations. First, we 
neglect momentum transfer from the dust to the gas caused 
by collisional drag. This allows us to decouple the equations 
of dust dynamics from the full system of model equations 
and to treat them independently. Second, we treat each dust 
species as an individual component neglecting grain destruc¬ 
tion. All the dust particles of a certain type within a cell of 
the hydrodynamic grid have the same velocity. Third, we 
assume that drifting dust grains affect neither the gas tem¬ 
perature nor the gas dynamics. To take into account dust 
absorption of UV radiation in the hydrodynamic and chemi¬ 
cal model of an H ll region, we use a ‘generalized’ dust which 
has a mean radius of 1.3 x 10“® cm and is mixed with gas 
with a uniform mass fraction of 1 per cent. Thus, our model 
is inconsistent in the dust treatment, but a much more com¬ 
plicated model would be needed to achieve consistency in 
this respect. 

The motion of charged grains in a magnetic field 
is also affected by the Lorenz force. In this paper, 
the interaction of charged grains with a magnetic field 
is not considered as this would require knowledge of 
the magnetic field geometry in the H ll region. The 
strength of th e magnetic field has been studie d obser - 
vationally by iHarvev-Smith. Madsen fc Gaensleil ll201lll . 
while its large-scale morphol ogy has been modelled by a 
numb e r of authors, including [Krumholz . Stone fc Gar diner! 
(12001^1: iG endele v k, Krumholj ( 2012l ~) ; [Arthur et al. ( 201lli ; 
iMackev &: LimI "( 20111 1. If the large-scale magnetic field is 
regular, we may expect that dust grains would migrate 
mostly along the magnetic field lines, which would lead 
to oblate or arc-like morphologies of dust distribution. So, 
the radial dust density profiles presented in the next sec¬ 
tion could be mostly relevant for regular magnetic field ge¬ 
ometries with some factor depending on the magnetic field 
strength and orientation. A magnetic field with a signifi¬ 
cant irregular component requires much more sophisticated 
modelling. 


Table 4. Adopted dust types. 



Gomposition 

Radius (A) 

Polycyclic aromatic hydrocarbons 

- 

5 

Very small grains (VSGs) 

graphite 

30 

Intermediate size grains (ISGs) 

graphite 

200 

Big grains (BGs) 

silicate 

3000 


3 DUST DRIFT IN AN EXPANDING H ii 

REGION 

In this study we do not attempt to describe in de¬ 
tail any particular object, but nevertheless some fidu¬ 
cial parameters are needed to model an expanding Hll 
region. Thus, we adopted the observed data for the 
Hii region of RC W 120, which has an almost per¬ 
fect ring-like shape ( Deharveng. Zjwa gno fc CaplarJ l2005l : 
iDeharveng et al.ll200£ - Izavagno et al.ll2010l l. so it is a good 
prototype for our simulation. Its parameters are given in 
Table 11 

In our model, the configuration that approximately cor¬ 
responds to RCW 120 is reached in about 630 kyr after the 
onset of ionization. The distributions of ionized, atomic and 
molecular hydrogen at this time are shown in Fig. [T] (top 
panel). 

The hydrogen number density within the Hll region 
is about 40 cm“®, half the observational estimate for 
RCW 120 (see Table d. The region is surrounded by a 
dense neutral envelope, which was swept-up by the expand¬ 
ing shock front. The enhanced density in the compressed 
envelope is about five times the density in the molecular re¬ 
gion. The thermal structure of the region is illustrated in 
the middle panel of Fig. [T] 

We present results for four grain sizes, three of which 
(PAHs, VSCs, and BCs) correspond to ‘conventional’ dust 
grains that are typically considered in various dust stud¬ 
ies. We also consider grains that are bigger than VSCs but 
smaller than BCs. These intermediate-sized grains (ISCs) 
are added to the model to estimate possible outcomes of a 
more detailed treatment of the dust size distribution. The 
parameters of the dust types studied are summarized in Ta¬ 
ble H 

The grain charge is determined by the balance between 
competing processes of electron and ion accretion and photo¬ 
electric emission. Illumination conditions inside the Hll re¬ 
gion mostly drive grains to positive charges (see the bottom 
panel of Fig. [T]). For reference, maximum grain charges for 
13.6 eV photon absorption, such that the next photoelectric 
emission is prevented by the electrostatic barrier, are -|-3, 
+19, +128 and +1167 e for PAHs, VSCs, ISCs and BCs, 
respectively (Eq. (22) from IWeingartner fc Draind [ 200 II I . 
The further away we move from the star, the less effective 
photoelectric emission becomes due to radiation absorption 
and dilution. At the same time, the electron accretion rate 
is pretty much the same over the entire Hll region. This 
is why the mean charge within the ionized region decreases 
outwardly and reaches negative values for PAHs, VSCs and 
BCs. ISCs are somewhat special in this respect as the neg¬ 
ative gradient of their charge is much shallower than for 
other grain types, and they do not become negative at all. 
The reason is that dust grains of this size absorb photons 
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Table 3. Parameters of RCW 120 used in the modelling. 


Quantity 

Value 

Reference 

Tstar 

Rstax 

35000 K 

6Rq 

3 X 10® cm“® 

Georeelin & Georeelin ('1970); Gramoton GQ?!): Martins et al. ('2010'): Diaz-Miller. Franco & Shore fl998) 
Diaz-Miller. Franco & Shore^^lOO^i 

n(H -I- H2)init 

Zavaeno et al.^^^OO^ 

n(H+) 

86 cm“® 

Zavaefno et al.^^^OO^ 

-^Hllregion 

4' or 1.5 pc 

Adopting distance of 1.3 kpc from^uss^ (200^ 




0 0.5 1 1.5 2 

Distance from the star [pc] 


Figure 1. Top panel: Number density of ionized and neutral gas 
components for t = 630 kyr. Middle panel: Gas and dust tem¬ 
peratures. Bottom panel: Mean charges of PAHs, very small and 
intermediate graphite grains, and big silicate grains. The bottom 
panel is split into two parts to present positively and negatively 
charged grains. 


(and emit electrons) very effectively. The negative gradient 
for grain charges is quite important for the resulting dust 
density profiles, as will be discussed later. The sharp charge 
rise at the edge of the ionized region is caused by a sharp 
decrease in the electron number density, and the dip located 
further away is related to the enhanced electron accretion 
rate at the density bump. 

As details of dust drift sensitively depend on a grain 
charge, below we present results for neutral and charged 
grains separately. Also we will check how important a spe¬ 
cific way of taking the grain charge into account is. Possible 
options are to consider the mean charge only for grains of a 
given size, the charge distribution for an ensemble of single¬ 
size grains or the fluctuating charge for a grain. These ap¬ 
proaches produce results with different levels of detail (and 
computational cost), so it may be desirable to find the sim¬ 
plest necessary approach. 


3.1 Neutral grains 

Grains within an ff ll region are definitely not neutral, so we 
present a calculation of dust drift for neutral PAHs, VSGs, 
ISGs and BGs only as a test case. It is assumed that initially 
all the dust components are well mixed with the gas. When 
the expansion of the Hll region starts, radiation pressure 
pushes grains away from the star, and the only force that 
counteracts their motion is the collisional drag. As the gas 
density is not very high, grain velocities for all the four grain 
types significantly exceed the gas velocity as shown in the 
bottom left-hand panel of Fig. (2] being greater than 100 km 
s“^ in the immediate star vicinity. The difference becomes 
smaller than 10 km s”*^ only near the border of the ionized 
region. 

The bigger the particle, the weaker it is coupled to the 
gas, so ISGs and BGs are swept out of the Hll region com¬ 
pletely (their density profiles coincide on the graph), while 
PAHs retain their initial abundance in the outer part of the 
region (remember that dust destruction is not taken into 
account), as seen in the top left-hand panel of Fig. [51 By 
630 kyr, about 25 per cent and 50 per cent of the Hll re¬ 
gion volume is free from PAHs and VSGs, respectively. The 
gas and dust velocities are the same in the compressed neu¬ 
tral envelope, where the dust is coupled to the gas. On this 
and other similar plots, all the dust densities are normalized 
to 1 per cent of the gas number density. This normaliza¬ 
tion emphasizes the difference between density profiles and 
drag regimes for grains of different sizes. Their actual rel¬ 
ative and absolute abundances should be estimated using 
emission modelling and other constraints set by general dust 
models. This will be presented in a subsequent study. 
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Neutral grains 
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Figure 2. Top panels: Density profiles of neutral (left) and charged (right) PAHs, VSGs, ISGs and BGs normalized to 1 per cent of the 
gas number density. Bottom panels: Radial velocities of neutral (left) and charged (right) PAHs, VSGs, ISGs and BGs along with gas 
velocity. The time since the beginning of the expansion is 630 kyr. 


3.2 Grains with mean charge 

The dynamics of charged grains differs drastically from that 
of neutral grains due to additional Coulomb drag. As a first 
step in the drift modelling of charged grains, we consid¬ 
ered the drift of grains with a representative mean charge. 
Whether the representative charge leads to a representative 
density profile will be discussed later. Velocities and spatial 
distribution of PAHs, VSGs, ISGs and BGs for t = 630 kyr 
are shown on the right-hand panels of Fig. [2] (their mean 
charges are plotted on the bottom panel of Fig. [T]). Obvi¬ 
ously, the spatial distributions of charged grains are not the 
same as those of neutral ones. Coulomb drag notably slows 
down the grain motion, making the dust more coupled to 
the gas. The velocities of charged PAHs, VSGs and ISGs 
only slightly exceed the gas velocity in most of the Hll re¬ 
gion. This is why charged particles, unlike neutral particles, 
stay relatively abundant in the ionized gas. On the other 
hand, dust drift due to radiation pressure has a profound 
effect on both charged and neutral BGs, significantly de¬ 
creasing their density in the ionized gas in both cases, even 
though the abundance of charged BGs exceeds significantly 
the abundance of BGs in the test case of zero charges. 

The velocities of all grains become greater than 100 km 
s“^ within about 0.1 pc from the star. Thus, the vicinity 
nearest to the star is cleaned from dust completely. Note 
the discontinuity in the velocity profile of BGs close to the 
star, which appears because the equation of motion with 
the Goulomb term may have three solutions for the veloc¬ 
ity Id raind lioilal ). In the general case, it has two stable 
solutions (low and high velocity) and one unstable solution 
(intermediate velocity). In most of the Hll region, we pick 
the stable low velocity solution as being more physical, but 


this solution disappears in the star’s vicinity, so the high 
velocity solution is presented. 

There are two drift regimes depending on whether there 
is a notable zone where the grains have zero or close to zero 
charge (PAHs and VSGs). In such regions, the Coulomb drag 
force vanishes and the grains fly outside faster till they en¬ 
ters the region where their charge becomes negative and 
where the Coulomb drag force again becomes significant. 
This leads to a pile-up of PAHs and VSGs at 0.8 and 1.3 pc, 
respectively. For bigger grains, the zero charge zone is ab¬ 
sent (ISGs) or has much smaller extent (BGs), so another 
drift regime operates. While the radiation pressure force de¬ 
creases outside, which slows the grain drift, the grain charge 
also decreases, lowering the grain coupling to the gas and 
increasing its velocity. For the considered physical parame¬ 
ters, Flp/Pdrag increases outside, so the grain velocity also 
increases outside. As an outer part of the ionized zone is 
cleared more effectively from grains compared to the inner 
part, a dust density bump emerges inside the Hll region. 

So, we see that when we take the mean charge into ac¬ 
count, a dust density bump emerges inside the Hll region. 
This is an intriguing finding in view of the double-ring ge¬ 
ometries of Hll regions in the mid-IR. Let us now check 
whether these features are preserved when we introduce a 
more detailed description of grain charging. 

3.3 Grains with charge distribution 

So far we have assumed that all the grains of a certain type 
have the same charge at a given time and at a given location. 
However, grains in the ensemble of equal-sized grains will 
attain slightly different charges, so that there will be a charge 
distribution with a non-zero dispersion. The dispersion is 
proportional to the charge value, so bigger grains have larger 
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charge scatter (see eqs. (5.81)-(5.82) in lTieler3 l2005h . We 
use two approaches to take into acconnt the impact of the 
charge dispersion on the dust density profiles. 

In the first approach, we consider an ensemble of single¬ 
type grains with a distribution of charges, which is charac¬ 
terized by a standard deviation for a given location and 
time. We consider 37 bins for grain charge, assuming that 
each bin contains the same fraction of the total dnst mass, 
and trace grain motion for each bin separately. In this ap¬ 
proach, the charge of a grain changes with time in accord 
with the mean charge and standard deviation for this grain 
type, so that Zi = Z + fia^ with factors fi 7^ fi{t) allow¬ 
ing an even mass distribution between bins. In the first ap¬ 
proach, grains are allowed to have a fractional charge as we 
try to trace the motion of the whole ensemble, rather than 
of individual grains. Density proHles for individual bins are 
shown in Fig. |3]as rainbow-like groups of lines. The sum of 
individual profiles is shown by the thick grey line. For com¬ 
parison, we also show the density proHles for grains with the 
mean charge (black lines). 

Density radial proHles for ISGs and BGs are nearly the 
same as when all the grains have the same (mean) charge. 
The depression in the ISGs density radial profile on the inner 
side of the dense shell becomes deeper for grains having a 
charge smaller than the mean. Similar behaviour is observed 
for BGs as well. The density in the ionized region is lower 
for grains with lower charge and is higher for grains with 
higher charge. But they all sum up to nearly the same total 
density as when there is no charge dispersion. 

On the other hand, the effect of the grain charge dis¬ 
tribution for PAHs and VSGs is profound. While density 
profiles for grains in each charge bin show the same general 
morphology, with the inner density bumps and the pile-up 
peaks, the location of the pile-up peak is a strong function of 
the grain charge. The peak is located much closer to the star 
for grains at the lower boundary of the considered charge 
range and is nearly absent for grains at the upper boundary 
of the range. The net result is that the actual density profile 
is quite flat, with all the pile-up peaks spread evenly over 
the ionized region. 

In the approach, we have just described, we assume that 
the grain charge changes only in response to changes in the 
physical environment in the sense that both Z and change 
with time, but fi stays constant. But the charge distribu¬ 
tion arises not only because there are grains with different 
charges, but also because the charge on each particular grain 
fluctuates with time. In the second approach, we consider 
grains with a fluctuating charge, which is set to a Gaussian 
random value with a given mean and standard deviation at 
every hydrodynamic time step (only integer charge values 
are allowed). Five realizations are considered with different 
random seeds. The results are shown in Fig. |4]for ISGs and 
BGs and in Fig. [5] for PAHs and VSGs. 

Again, the results for ISGs and BGs are nearly the 
same as when all the grains are assumed to have the mean 
charge determined by the current physical environment. 
Their mean charge is large, and its fluctuations are not sig¬ 
nificant, so that the dust density profiles for the three con¬ 
sidered cases do not differ much. 

The situation is somewhat different for VSGs. As we 
mentioned above, a realistic grain charge distribution tends 
to smooth the resulting density profile in comparison with 


when all the VSGs are assumed to have the mean charge. 
This tendency is preserved also when each grain is assumed 
to have a fluctuating charge. It does not really matter how 
exactly the charge distribution is taken into account for 
VSGs. Even though charge excursions from the mean for 
a single grain can be quite significant (bottom right-hand 
panel of Fig. El. the resulting density profile (top right-hand 
panel of Fig. [5] red band) is very similar to the profile for 
the grain ensemble with a non-fluctuating charge distribu¬ 
tion (top right-hand panel of Fig. O grey band). 

Finally, PAHs behave differently in all three considered 
cases. Their mean charge can attain both positive and neg¬ 
ative values, ensuring effective coupling between PAHs and 
the gas over most of the H ll region and causing the pile-up 
peak to appear near the distance range where the charge is 
zero. In the model with the grain ensemble, positively and 
negatively charged PAHs are mixed with each other, and all 
the density profile features are smeared out. The situation 
changes dramatically when we add charge fluctuations for 
each particular PAH particle. Their mean charge is quite 
small, so the charge of a PAH particle fluctuates between 
-2 and -|-3, being neutral for a significant portion of the 
time (bottom left-hand panel of Fig. [5|. As a result, the dis¬ 
tance range in which the PAHs are decoupled from the gas 
is significantly wider than when all the PAHs are assumed 
to have the mean charge. An inner cavity is formed due to 
the effective escape of PAHs from the star’s vicinity when 
PAH particles have zero charge. 

We conclude that the mean charge is representative in 
terms of dust dynamics for larger grains (ISGs and BGs), 
but a more sophisticated consideration is needed for smaller 
grains (PAHs and VSGs). For PAHs, a fluctuating grain 
charge seems to be the only realistic approach among the 
approaches we have considered. This approach is preferable 
for VSGs also. 

3.4 Charged grains: chemical composition 

The impact of grain chemical composition on the dust drift 
is twofold. First, chemical composition affects grain charge 
both through different absorbing efficiencies and photoelec¬ 
tron work functions for graphites and silicates. Generally, as 
graphites have larger absorbing coefficients and lower pho¬ 
toelectron work function, they have larger charges and are 
more efficiently coupled to the gas. Second, this coupling is 
compensated somewhat by the stronger radiation pressure 
acting on the more absorbing graphitic grains. The overall 
effect of grain chemical composition on dust drift is shown 
in Fig.[^for the case when all the grains are assumed to have 
fluctuating charges. Graphitic grains tend to have more ex¬ 
tended distributions and are better coupled to the gas in 
comparison with silicate grains, though the difference is not 
qualitative. 


4 DISCUSSION 

The presented study is our next _ step (following 

IPavlvuchenkov. Kirsanova fc Wiebel l2013l l towards a 
consistent interpretation and explanation of observed IR 
intensity distributions in Hll regions around young stars. 
There is a general consensus that 8 /rm emission is formed 
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Figure 3. Density profiles of charged PAHs, VSGs, ISGs and BGs for t = 630 kyr normalized to 1 per cent of the gas number density. 
Golour lines illustrate density variations due to charge dispersion within an ensemble of single-sized grains. 
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Figure 4. Comparison of approaches for grain charge evolution. Top panels: Density profiles for ISGs and BGs normalized to 1 per cent 
of the gas number density. Red lines show density profiles for grains with a fluctuating charge. The thick grey line and black lines are 
the same as in Fig. 0 Bottom panel: ISGs and BGs charges (dots between upper and lower sub-panels represent zero charges). 
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Figure 5. Same as Fig. [4| but for PAHs and VSGs. Fluctuations in the density profile of PAHs (red lines) reflect dispersion due to 
different random seeds in the stochastic charge modelling. For other grain types, these fluctuations are smaller than the line width. 
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Figure 6. Top panels: Density profiles of graphitic and silicate VSGs, ISGs and BGs normalized to 1 per cent of the gas number density. 
Bottom panels: Their corresponding mean charges. 


by PAHs, so that the ring-like 8 fim emission maps hint at 
the absence of PAHs in the interiors of Hll regions, where 
they are either destroyed or blown away by the stellar wind 
or by radiation pressure. 

The presence of 24 fim emission in the interiors of H ll 


regions indicates that grains larger than PAHs may exist 
there, but the exact origin of the 24 ^tm e mission is still a 
puzzle. According to IPaladini et all ||2012| '). at A > 70 /rm 
IR spectra of 16 Hll regions up to 100 /rm, obtained with 
Herschel, can be generated by big grains with Tdust about 
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20-30 K, while smaller grains with Tdust about 50-90 K are 
responsible for IR emission with A < 70 ^m. 

Dynamical models of H ll regions s how the impor tance 
of stellar win ds, see e.g. reviews by lArthurl (l2007l l and 
iHennevI l|2007l ~l. However, it is not easy to explain how small 
grains may survive in a win d-blown bubble. This has led 
lEverett fc Churchwelll (l2010l l to the conclusion that dust 
within an IR bubble may be constantly replenished from 
evaporating protopla netary discs (dusty clou dlets) embed¬ 
ded in the H ll region. lOchsendorf et al.l (l2014h related 24 pm 
emission to the material that is photo-evaporated from the 
inner bubble walls and then is shaped into a structure, which 
they coined a dust wave, by the radiation pressure of an 
exciting star. An important feature of the latter model is 
the pressure leakage from the bubble into the surrounding 
medium. 

We explore whether dust distribution within an IR bub¬ 
ble can be explained without additional assumptions, like 
disc or inner wall evaporation, but simply due to dust evo¬ 
lution within a pressure-driven bubble. Specifically, in this 
paper we consider the effect of radiatively driven dust drift 
on the spatial distribution of different dust types. Our con¬ 
clusion is that the combined action of radiation pressure and 
Coulomb drag on charged dust grains results in a signihcant 
stratification of dust of different sizes and thus can be a po¬ 
tential explanation for the observed IR morphology of H ll 
regions. This needs to be further supported by radiation 
transfer simulations, which will be presented in forthcoming 
papers. 

Also, we have to keep in mind that there are other phys¬ 
ical processes that are relevant in explaining the appearance 
of IR bubbles, for example, radiation pressure on t he gas. 
The sta tionary models of an Hll region presented bv lDraina 
|201 i 3), which account for the radiation pressure on the gas, 
demonstrate that density profiles inside the H ll region may 
have very different shapes depending on the model param¬ 
eters. In particular, there are solutions with inner cavities 
that may be responsible for the dehcit of dust in the inte¬ 
riors of Hll regions. The importance of radiation pressure 
on the gas and/or stellar wind is supported by the spa¬ 
tial correlation between 20 cm free-free emission and 24 pm 
dust emission towards s ome IR bubbles (e.g., N36 and N49; 
iDeharveng et al.l (l201Cll H. 

Neither radiation pressure on the gas nor the stellar 
wind is taken into account in the presented model, and the 
density distribution of the ionized gas is nearly flat. Both 
factors, acting in concert with the radiation pressure on the 
dust, may be responsible for the absence of dust within the 
24 pm ring. We plan to include the radiation pressure on the 
gas with the full dynamical coupling between the gas and 
dust it our future studies. But for the moment, our results 
indicate that a non-uniform dust distribution can appear 
even without these factors. 

Another effect that may be important is the recy¬ 
cling/evolution of dust triggered by UV radi ation a n d hig h 
temperature in Hii regions. According to I Jones! ll2009f) . 
the aromatisation and subsequent fragmentation of hydro¬ 
genated amorphous carbon grains can be a primary source 
of the aromatic emission band carriers. With this in mind, 
one may suggest that PAHs do not even need to be ini¬ 
tially present in the molecular gas. The 8 pm emission can 
be generated by the transient layer of dust particles that 


are formed and quickly destroyed just behind the ionization 
front, rather than by PAHs that have been somehow re¬ 
moved from the ionized region. To check this idea, it is nec¬ 
essary to develop a dynamical model coupled with a model 
of dust microevolution. 


5 CONCLUSIONS 

In this paper we investigated dust dynamics in a pressure- 
driven expanding H ll region to understand the appearance 
of IR bubbles. As drag force and radiation pressure force 
strongly depend on dust size, their combined action leads 
to dust differentiation within an Hll region. We particu¬ 
larly focused on the effect of radiation pressure on charged 
grain motion and do not consider stellar wind and dust de¬ 
struction. The grain charge is consistently evaluated to the 
compute aerodynamic drag force. Our conclusions can be 
summarized as follows: 

• Grain charge is of utter importance for dust velocity 
and t he resultant density proHles (see also lGail fc Sedlmavil 
I 1979 IH . If all the grains are assumed to be neutral, the inner 
parts of the H ll region are almost devoid of dust. The grain 
charge increases dust friction with the gas and suppresses 
dust removal from the H ll region. 

• The grain ensemble is characterized by the mean charge 
and charge dispersion. It is crucial to consider both param¬ 
eters in the dynamics of small grains (a < 100 A). Goulomb 
drag for larger grains (a > 100 A) can be computed using 
their mean charge. 

• The interplay between radiation pressure and drag gen¬ 
erally leads to double-peak dust density prohles. The outer 
peak corresponds to the gas density maximum on the ion¬ 
ization and dissociation fronts. The inner peak arises due 
to the negative gradient of grain charge within the ionized 
region. Grains in the outer parts of the region have smaller 
charges and are less coupled to the gas. Sweeping of these 
grains leads to a density dip between the inner ionized region 
and the outer density maximum. The details of this process 
strongly depend on grain size. 

• PAHs, having the largest surface-to-mass ratio, are well 
coupled to the gas whenever they have a non-zero charge. 
PAHs with zero charge are easily blown away from the star’s 
vicinity. Near-zero charge fluctuations during the expansion 
lead to removal of PAHs from the inner part of the Hll 
region. So, in principle, the lack of 8 pm emission inside the 
Hii region can be explained not only by the destruction of 
PAHs but also by the radiative removal of PAHs if their 
charge fluctuates around zero. 

• Very small grains acquire zero charge less frequently 
than PAHs and, consequently, are better coupled to the gas 
and do not escape ionized region easily. The density distri¬ 
bution of VSGs is the most uniform relative to other grain 
types. Stochastically heated VSGs may partially (along with 
ISGs) explain the double-peak morphology of IR bubbles at 
24 pm. 

• The intriguing feature of intermediate-sized grains {a ~ 
200 A) is that their emission could contribute to both mid- 
and far-IR bands. The inner peak on the ISG density distri¬ 
bution is several orders of magnitude smaller than the outer 
peak, but the negative dust temperature gradient may lead 
to an inverse ratio of corresponding intensities at 24 pm. At 
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the same time, the outer peak in the ISG distribution may 
be observable at 70-100 pm. 

• Big grains are subject to the strongest sweeping by the 
radiation pressure. Being able to contribute to far-IR, they 
are worthy candidates for the outer ring for emission at 70- 
100 pm. 

The real dust size distribution in H ll regions is more com¬ 
plicated than described by the four dust types listed above. 
However, consideration of PAHs, VSGs, ISGs and BGs is 
necessary to understand the contribution of different sub¬ 
samples of the grain size distribution in the dust emission 
spectrum. While such problems are generally ill defined, ad¬ 
ditional information can be retrieved from the Hll region 
morphology accompanied by grain dynamics modelling. 
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APPENDIX A: GRAIN GHARGE 


The grain charge is primarily determined by photoelectric 
emission and accretion of ions and electrons. The proba¬ 
bility f{Zgr) of finding a grain having charge ZgrC can be 
deter mined from the equation of detailed balance (iTielensI 

I2OO5II : 


/(^gr) [>^pe(^gr) + Jion(^gr)] — 

~ f{Zgr + i)JeiZgr + 1), (Al) 

where Jpe [s~^] is a photoelectric emission rate, Je and 
Tion are the electron and ion accretion rates, respec¬ 

tively. Writing down the above equation for every •^gri one 
finds an equation for /(Zgr): 

fiZg.) = fo n (A2) 

'T- I 1 Je[Z) 

z=Zq + 1 ^ ' 

for Zgr > Zq and 


Zo-i 

fiZgr) = fo n 
z = Zgr 


Je{z + 1) 

</ion(-2) 


(A3) 


for Zgr < Zo- Here /o = fiZo). From the computational 
point of view, it is important that Zo should correspond to 
the maximum of f{Zgr), which can be estimated from the 
equation: 


Jpe(^o) ~ MZo). (A4) 

Equations (IA2ll - (IA3ll are closed by 


E = (A5) 


Photon attenuation depth and electron mean free path are 
assumed to be equal to lOOAand lOA, respectively. 

The accretion rate of part icles with number den sity n 
and mass m was adopted from lPraine fc SutinI (Il987l l: 


J acc — ri 


8kT, 


gas 2 
TT Ugj- 


J{c 


kr 1 ^gas , Zgi ) 


(A6) 


Here J accounts for Goulomb focus ing. As sticking proba¬ 
bilitie s are not well known, we follow IWeingartner fc Draind 
(l200lll in assuming sticking probabilities for the electrons 
and ions to be equal to 0.5 and 1, respectively. 

The photoemission rate depends on the mean intensity 
Ji, [erg/(cm^ s Hz sr) ] and photo-ionization yield Y: 


Jpe — 


47rJ,, 

hv 


Gabs Ugr )T (zz, ^gr , Ugr ) dv 


(A7) 


and is calculated according to lWeingartner fc Draind l|200j l. 


REFERENCES 

Aannestad P. A., Emery R. J., 2001, A&A, 376, 1040 
Aitken D. K., Griffiths J., Jones B., 1977, MNRAS, 179, 
179 

Anderson L. D. et al., 2012, A&A, 542, AlO 
Arthur S. J., 2007, Wind-Blown Bubbles around Evolved 
Stars, Hartquist T. W., Pittard J. M., Falle S. A. E. G., 
eds., p. 183 


© 2015 RAS, MNRAS OOO.fTWTTI 






















12 V.V. Akimkin, M.S. Kirsanova, Ya.N. Pavlyuchenkov, D.S. Wiehe 


Arthur S. J., Henney W. J., Mellema G., de Colie F., 
Vazquez-Semadeni E., 2011, MNRAS, 414, 1747 
Arthur S. J., Kurtz S. E., Franco J., Albarran M. Y., 2004, 
ApJ, 608, 282 

Chini R., Kruegel E., Wargau W., 1987, A&A, 181, 378 
Churchwell E. et ah, 2006, ApJ, 649, 759 
Compiegne M. et ah, 2011, A&A, 525, A103 
Crampton D., 1971, AJ, 76, 260 
Deharveng L. et ah, 2010, A&A, 523, A6 
Deharveng L., Zavagno A., Caplan J., 2005, A&A, 433, 565 
Deharveng L., Zavagno A., Schuller F., Caplan J., Pomares 
M., De Breuck C., 2009, A&A, 496, 177 
Diaz-Miller R. L, Franco J., Shore S. N., 1998, ApJ, 501, 
192 

Draine B. T., 2011a, ApJ, 732, 100 

Draine B. T., 2011b, Physics of the Interstellar and Inter- 
galactic Medium. Princeton University Press, 2011. ISBN: 
978-0-691-12214-4 

Draine B. T., Li A., 2007, ApJ, 657, 810 
Draine B. T., Salpeter E. E., 1979, ApJ, 231, 77 
Draine B. T., Sutin B., 1987, ApJ, 320, 803 
Everett J. E., Churchwell E., 2010, ApJ, 713, 592 
Ferland G. J. et ah, 2013, Rev. Mex. Ast., 49, 137 
Gail H. P., Sedlmayr E., 1979, A&A, 77, 165 
Gendelev L., Krumholz M. R., 2012, ApJ, 745, 158 
Georgelin Y. P., Georgelin Y. M., 1970, A&AS, 3, 1 
Gillett F. C., Forrest W. J., Merrill K. M., Soifer B. T., 
Capps R. W., 1975, ApJ, 200, 609 
Harper D. A., Low F. J., 1971, ApJL, 165, L9 
Harvey-Smith L., Madsen G. J., Gaensler B. M., 2011, ApJ, 
736, 83 

Henney W. J., 2007, How to Move Ionized Gas: An Intro¬ 
duction to the Dynamics of HII Regions, Hartquist T. W., 
Pittard J. M., Falle S. A. E. G., eds., p. 103 
Hollenbach D., McKee C. F., 1979, ApJS, 41, 555 
Hollenbach D., McKee C. F., 1989, ApJ, 342, 306 
Hosokawa T., Inutsuka S.-i., 2006, ApJ, 646, 240 
Jones A. P., 2009, in Astronomical Society of the Pacific 
Conference Series, Vol. 414, Cosmic Dust - Near and Far, 
Henning T., Griin E., Steinacker J., eds., p. 473 
Jonkheid B., Faas F. G. A., van Zadelhoff G.-J., van 
Dishoeck E. F., 2004, A&A, 428, 511 
Kirsanova M. S., Wiebe D. S., Sobolev A. M., 2009, As¬ 
tronomy Reports, 53, 611 

Krumholz M. R., Matzner C. D., 2009, ApJ, 703, 1352 
Krumholz M. R., Stone J. M., Gardiner T. A., 2007, ApJ, 
671, 518 

Kurucz R. L., 1979, ApJS, 40, 1 
Mackey J., Lim A. J., 2011, MNRAS, 412, 2079 
Martins F., Pomares M., Deharveng L., Zavagno A., Bouret 
J. C., 2010, A&A, 510, A32 
Mathews W. G., 1965, ApJ, 142, 1120 
Mathews W. G., 1967, ApJ, 147, 965 
Millar T. J., Farquhar P. R. A., Willacy K., 1997, A&AS, 
121, 139 

Nakano M., Kogure T., Sasaki T., Mizuno S., Sakka K., 
Wiramihardja S. D., 1983, Ap&SS, 89, 407 
Neufeld D. A., Kaufman M. J., 1993, ApJ, 418, 263 
Ney E. P., Allen D. A., 1969, ApJL, 155, L193 
Ochsendorf B. B., Verdolini S., Cox N. L. J., Berne O., 
Kaper L., Tielens A. G. G. M., 2014, A&A, 566, A75 
O’Dell C. R., Hubbard W. B., 1965, ApJ, 142, 591 


O’Dell C. R., Hubbard W. B., Peimbert M., 1966, ApJ, 
143, 743 

Paladini R. et ah, 2012, ApJ, 760, 149 
Panagia N., 1974, ApJ, 192, 221 

Pavlyuchenkov Y. N., Kirsanova M. S., Wiebe D. S., 2013, 
Astronomy Reports, 57, 573 
Petrosian V., Silk J., Field G. B., 1972, ApJL, 177, L69 
Russeil D., 2003, A&A, 397, 133 
Simpson R. J. et ah, 2012, MNRAS, 424, 2442 
Spitzer L., 1978, Physical processes in the interstellar 
medium. New York Wiley-Interscience, 1978. 333 p. 

Stein W. A., Gillett F. C., 1969, ApJL, 155, L197 
Stone J. M., Norman M. L., 1992, ApJS, 80, 753 
Tenorio-Tagle G., 1974, Ap&SS, 26, 111 
Tielens A. G. G. M., 2005, The Physics and Chemistry 
of the Interstellar Medium, Cambridge, UK: Cambridge 
University Press. Cambridge University Press 
Tielens A. G. G. M., Hollenbach D., 1985, ApJ, 291, 722 
van Dishoeck E. F., 1988, in Astrophysics and Space Sci¬ 
ence Library, Vol. 146, Rate Coefficients in Astrochem- 
istry, Millar T. J., Williams D. A., eds., pp. 49-72 
van Dishoeck E. F., Jonkheid B., van Hemert M. C., 2006, 
Faraday Discussions, 133, 231 
Watson C. et ah, 2008, ApJ, 681, 1341 
Weingartner J. C., Draine B. T., 2001, ApJS, 134, 263 
Wright E. L., 1973, ApJ, 185, 569 

Zavagno A., Pomares M., Deharveng L., Hosokawa T., Rus¬ 
seil D., Caplan J., 2007, A&A, 472, 835 
Zavagno A. et ah, 2010, A&A, 518, L81 


© 2015 RAS, MNRAS OOP.IHITTI 


